Evolutionary conserved microRNAs are ubiquitously expressed compared to tick-specific miRNAs in the cattle tick Rhipicephalus (Boophilus) microplus

Background: MicroRNAs (miRNAs) are small non-coding RNAs that act as regulators of gene expression in eukaryotes modulating a large diversity of biological processes. The discovery of miRNAs has provided new opportunities to understand the biology of a number of species. The cattle tick, Rhipicephalus (Boophilus) microplus, causes significant economic losses in cattle production worldwide and this drives us to further understand their biology so that effective control measures can be developed. To be able to provide new insights into the biology of cattle ticks and to expand the repertoire of tick miRNAs we utilized Illumina technology to sequence the small RNA transcriptomes derived from various life stages and selected organs of R. microplus. Results: To discover and profile cattle tick miRNAs we employed two complementary approaches, one aiming to find evolutionary conserved miRNAs and another focused on the discovery of novel cattle-tick specific miRNAs. We found 51 evolutionary conserved R. microplus miRNA loci, with 36 of these previously found in the tick Ixodes scapularis. The majority of the R. microplus miRNAs are perfectly conserved throughout evolution with 11, 5 and 15 of these conserved since the Nephrozoan (640 MYA), Protostomian (620MYA) and Arthropoda (540 MYA) ancestor, respectively. We then employed a de novo computational screening for novel tick miRNAs using the draft genome of I. scapularis and genomic contigs of R. microplus as templates. This identified 36 novel R. microplus miRNA loci of which 12 were conserved in I. scapularis. Overall we found 87 R. microplus miRNA loci, of these 15 showed the expression of both miRNA and miRNA* sequences. R. microplus miRNAs showed a variety of expression profiles, with the evolutionary-conserved miRNAs mainly expressed in all life stages at various levels, while the expression of novel tick-specific miRNAs was mostly limited to particular life stages and/or tick organs. Conclusions: Anciently acquired miRNAs in the R. microplus lineage not only tend to accumulate the least amount of nucleotide substitutions as compared to those recently acquired miRNAs, but also show ubiquitous expression profiles through out tick life stages and organs contrasting with the restricted expression profiles of novel tick-specific miRNAs

Evidence of a tick RNAi pathway by comparative genomics and reverse genetics screen of targets with known loss-of-function phenotypes in Drosophila

Background The Arthropods are a diverse group of organisms including Chelicerata (ticks, mites, spiders), Crustacea (crabs, shrimps), and Insecta (flies, mosquitoes, beetles, silkworm). The cattle tick, Rhipicephalus (Boophilus) microplus, is an economically significant ectoparasite of cattle affecting cattle industries world wide. With the availability of sequence reads from the first Chelicerate genome project (the Ixodes scapularis tick) and extensive R. microplus ESTs, we investigated evidence for putative RNAi proteins and studied RNA interference in tick cell cultures and adult female ticks targeting Drosophila homologues with known cell viability phenotype. Results We screened 13,643 R. microplus ESTs and I. scapularis genome reads to identify RNAi related proteins in ticks. Our analysis identified 31 RNAi proteins including a putative tick Dicer, RISC associated (Ago-2 and FMRp), RNA dependent RNA polymerase (EGO-1) and 23 homologues implicated in dsRNA uptake and processing. We selected 10 R. microplus ESTs with >80% similarity to D. melanogaster proteins associated with cell viability for RNAi functional screens in both BME26 R. microplus embryonic cells and female ticks in vivo. Only genes associated with proteasomes had an effect on cell viability in vitro. In vivo RNAi showed that 9 genes had significant effects either causing lethality or impairing egg laying. Conclusion We have identified key RNAi-related proteins in ticks and along with our loss-of-function studies support a functional RNAi pathway in R. microplus. Our preliminary studies indicate that tick RNAi pathways may differ from that of other Arthropods such as insects

Genome sequencing and analysis of the paclitaxelproducing endophytic fungus \u3cem\u3ePenicillium aurantiogriseum\u3c/em\u3e NRRL 62431

Background Paclitaxel (Taxol™) is an important anticancer drug with a unique mode of action. The biosynthesis of paclitaxel had been considered restricted to the Taxus species until it was discovered in Taxomyces andreanae, an endophytic fungus of T. brevifolia. Subsequently, paclitaxel was found in hazel (Corylus avellana L.) and in several other endophytic fungi. The distribution of paclitaxel in plants and endophytic fungi and the reported sequence homology of key genes in paclitaxel biosynthesis between plant and fungi species raises the question about whether the origin of this pathway in these two physically associated groups could have been facilitated by horizontal gene transfer. Results The ability of the endophytic fungus of hazel Penicillium aurantiogriseum NRRL 62431 to independently synthesize paclitaxel was established by liquid chromatography-mass spectrometry and proton nuclear magnetic resonance. The genome of Penicillium aurantiogriseum NRRL 62431 was sequenced and gene candidates that may be involved in paclitaxel biosynthesis were identified by comparison with the 13 known paclitaxel biosynthetic genes in Taxus. We found that paclitaxel biosynthetic gene candidates in P. aurantiogriseum NRRL 62431 have evolved independently and that horizontal gene transfer between this endophytic fungus and its plant host is unlikely. Conclusions Our findings shed new light on how paclitaxel-producing endophytic fungi synthesize paclitaxel, and will facilitate metabolic engineering for the industrial production of paclitaxel from fungi

Ezrin interacts with the SARS coronavirus spike protein and restrains infection at the entry stage

© 2012 Millet et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Background: Entry of Severe Acute Respiratory Syndrome coronavirus (SARS-CoV) and its envelope fusion with host cell membrane are controlled by a series of complex molecular mechanisms, largely dependent on the viral envelope glycoprotein Spike (S). There are still many unknowns on the implication of cellular factors that regulate the entry process. Methodology/Principal Findings: We performed a yeast two-hybrid screen using as bait the carboxy-terminal endodomain of S, which faces the cytosol during and after opening of the fusion pore at early stages of the virus life cycle. Here we show that the ezrin membrane-actin linker interacts with S endodomain through the F1 lobe of its FERM domain and that both the eight carboxy-terminal amino-acids and a membrane-proximal cysteine cluster of S endodomain are important for this interaction in vitro. Interestingly, we found that ezrin is present at the site of entry of S-pseudotyped lentiviral particles in Vero E6 cells. Targeting ezrin function by small interfering RNA increased S-mediated entry of pseudotyped particles in epithelial cells. Furthermore, deletion of the eight carboxy-terminal amino acids of S enhanced S-pseudotyped particles infection. Expression of the ezrin dominant negative FERM domain enhanced cell susceptibility to infection by SARS-CoV and S pseudotyped particles and potentiated S-dependent membrane fusion. Conclusions/Significance: Ezrin interacts with SARS-CoV S endodomain and limits virus entry and fusion. Our data present a novel mechanism involving a cellular factor in the regulation of S-dependent early events of infection.This work was supported by the Research Grant Council of Hong Kong (RGC#760208)and the RESPARI project of the International Network of Pasteur Institutes

Temperature alters gene expression in mosquitoes during arbovirus infection

ABSTRACTArthropod-borne viruses (arboviruses) such as dengue, Zika and chikungunya constitute a significant proportion of the global disease burden. The principal vector of these pathogens is the mosquito Aedes (Ae.) aegypti, and its ability to transmit virus to a human host is influenced by environmental factors such as temperature. However, exactly how ambient temperature influences virus replication within mosquitoes remains poorly elucidated, particularly at the molecular level. Here, we use chikungunya virus (CHIKV) as a model to understand how the host mosquito transcriptome responds to arbovirus infection under different ambient temperatures. We exposed CHIKV-infected mosquitoes to 18 °C, 28 °C and 32 °C, and found higher temperature correlated with higher virus replication levels, particularly at early time points post-infection. Lower ambient temperatures resulted in reduced virus replication levels. Using RNAseq, we found that temperature significantly altered gene expression levels in mosquitoes, particularly components of the immune response. The highest number of significantly differentially expressed genes in response to CHIKV was observed at 28 °C, with a markedly more muted effect observed at either lower (18 °C) or higher (32 °C) temperatures. At the higher temperature, the expression of many classical immune genes, including Dicer-2 in the RNAi pathway, was not substantially altered in response to CHIKV. Upregulation of Toll, IMD and JAK-STAT pathways was only observed at 28 °C. Time post infection also led to substantially different gene expression profiles, and this effect varied depending upon the which temperature mosquitoes were exposed to. Taken together, our data indicate temperature significantly modulates mosquito gene expression in response to infection, potentially leading to impairment of immune defences at higher ambient temperatures.</jats:p

Curated genome annotation of Oryza sativa ssp. japonica and comparative genome analysis with Arabidopsis thaliana

We present here the annotation of the complete genome of rice Oryza sativa L. ssp. japonica cultivar Nipponbare. All functional annotations for proteins and non-protein-coding RNA (npRNA) candidates were manually curated. Functions were identified or inferred in 19,969 (70%) of the proteins, and 131 possible npRNAs (including 58 antisense transcripts) were found. Almost 5000 annotated protein-coding genes were found to be disrupted in insertional mutant lines, which will accelerate future experimental validation of the annotations. The rice loci were determined by using cDNA sequences obtained from rice and other representative cereals. Our conservative estimate based on these loci and an extrapolation suggested that the gene number of rice is ~32,000, which is smaller than previous estimates. We conducted comparative analyses between rice and Arabidopsis thaliana and found that both genomes possessed several lineage-specific genes, which might account for the observed differences between these species, while they had similar sets of predicted functional domains among the protein sequences. A system to control translational efficiency seems to be conserved across large evolutionary distances. Moreover, the evolutionary process of protein-coding genes was examined. Our results suggest that natural selection may have played a role for duplicated genes in both species, so that duplication was suppressed or favored in a manner that depended on the function of a gene

Temperature modulates immune gene expression in mosquitoes during arbovirus infection.

The principal vector of dengue, Zika and chikungunya viruses is the mosquito Aedes aegypti, with its ability to transmit pathogens influenced by ambient temperature. We use chikungunya virus (CHIKV) to understand how the mosquito transcriptome responds to arbovirus infection at different ambient temperatures. We exposed CHIKV-infected mosquitoes to 18, 28 and 32°C, and found that higher temperature correlated with higher virus levels, particularly at 3 days post infection, but lower temperature resulted in reduced virus levels. RNAseq analysis indicated significantly altered gene expression levels in CHIKV infection. The highest number of significantly differentially expressed genes was observed at 28°C, with a more muted effect at the other temperatures. At the higher temperature, the expression of many classical immune genes, including Dicer-2, was not substantially altered in response to CHIKV. The upregulation of Toll, IMD and JAK-STAT pathways was only observed at 28°C. Functional annotations suggested that genes in immune response and metabolic pathways related to energy supply and DNA replication were involved in temperature-dependent changes. Time post infection also led to substantially different gene expression profiles, and this varied with temperature. In conclusion, temperature significantly modulates mosquito gene expression in response to infection, potentially leading to impairment of immune defences at higher temperatures

High-Throughput SuperSAGE for Digital Gene Expression Analysis of Multiple Samples Using Next Generation Sequencing

We established a protocol of the SuperSAGE technology combined with next-generation sequencing, coined “High-Throughput (HT-) SuperSAGE”. SuperSAGE is a method of digital gene expression profiling that allows isolation of 26-bp tag fragments from expressed transcripts. In the present protocol, index (barcode) sequences are employed to discriminate tags from different samples. Such barcodes allow researchers to analyze digital tags from transcriptomes of many samples in a single sequencing run by simply pooling the libraries. Here, we demonstrated that HT-SuperSAGE provided highly sensitive, reproducible and accurate digital gene expression data. By increasing throughput for analysis in HT-SuperSAGE, various applications are foreseen and several examples are provided in the present study, including analyses of laser-microdissected cells, biological replicates and tag extraction using different anchoring enzymes

De novo assembly of Euphorbia fischeriana root transcriptome identifies prostratin pathway related genes

Background Euphorbia fischeriana is an important medicinal plant found in Northeast China. The plant roots contain many medicinal compounds including 12-deoxyphorbol-13-acetate, commonly known as prostratin that is a phorbol ester from the tigliane diterpene series. Prostratin is a protein kinase C activator and is effective in the treatment of Human Immunodeficiency Virus (HIV) by acting as a latent HIV activator. Latent HIV is currently the biggest limitation for viral eradication. The aim of this study was to sequence, assemble and annotate the E. fischeriana transcriptome to better understand the potential biochemical pathways leading to the synthesis of prostratin and other related diterpene compounds. Results In this study we conducted a high throughput RNA-seq approach to sequence the root transcriptome of E. fischeriana. We assembled 18,180 transcripts, of these the majority encoded protein-coding genes and only 17 transcripts corresponded to known RNA genes. Interestingly, we identified 5,956 protein-coding transcripts with high similarity (>=75%) to Ricinus communis, a close relative to E. fischeriana. We also evaluated the conservation of E. fischeriana genes against EST datasets from the Euphorbeacea family, which included R. communis, Hevea brasiliensis and Euphorbia esula. We identified a core set of 1,145 gene clusters conserved in all four species and 1,487 E. fischeriana paralogous genes. Furthermore, we screened E. fischeriana transcripts against an in-house reference database for genes implicated in the biosynthesis of upstream precursors to prostratin. This identified 24 and 9 candidate transcripts involved in the terpenoid and diterpenoid biosyntehsis pathways, respectively. The majority of the candidate genes in these pathways presented relatively low expression levels except for 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase (HDS) and isopentenyl diphosphate/dimethylallyl diphosphate synthase (IDS), which are required for multiple downstream pathways including synthesis of casbene, a proposed precursor to prostratin. Conclusion The resources generated in this study provide new insights into the upstream pathways to the synthesis of prostratin and will likely facilitate functional studies aiming to produce larger quantities of this compound for HIV research and/or treatment of patients

Large-scale identification and characterization of alternative splicing variants of human gene transcripts using 56 419 completely sequenced and manually annotated full-length cDNAs

We report the first genome-wide identification and characterization of alternative splicing in human gene transcripts based on analysis of the full-length cDNAs. Applying both manual and computational analyses for 56 419 completely sequenced and precisely annotated full-length cDNAs selected for the H-Invitational human transcriptome annotation meetings, we identified 6877 alternative splicing genes with 18 297 different alternative splicing variants. A total of 37 670 exons were involved in these alternative splicing events. The encoded protein sequences were affected in 6005 of the 6877 genes. Notably, alternative splicing affected protein motifs in 3015 genes, subcellular localizations in 2982 genes and transmembrane domains in 1348 genes. We also identified interesting patterns of alternative splicing, in which two distinct genes seemed to be bridged, nested or having overlapping protein coding sequences (CDSs) of different reading frames (multiple CDS). In these cases, completely unrelated proteins are encoded by a single locus. Genome-wide annotations of alternative splicing, relying on full-length cDNAs, should lay firm groundwork for exploring in detail the diversification of protein function, which is mediated by the fast expanding universe of alternative splicing variants